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Abstract 

We present measurements of the branching fraction and time-dependent CP- violating 
(CPV) asymmetries in K^tt^ decays based on 227 milUon T{4S) — > BB decays 

collected with the BABAR detector at the PEP-II asymmetric-energy B factory at 
SLAC. We obtain a branching fraction K°tt°) = (11.4 ± 0.9 ± 0.6) • 10"^, 

the magnitude of the direct CPV asymmetry C^o^o = 0.06 it 0.18 ± 0.06 and the 
magnitude of the CPV asymmetry in the interference between mixing and decay 
^k°7t0 ~ ^-"^^-o ss i0.04, where the first error is statistical and the second systematic. 
All results are preliminary. 
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1 INTRODUCTION 



In the Standard Model (SM) the time-dependent CP symmetry violation in decays via 
h sqq {q G {u, d, s, c}) transitions is governed by the Cabibbo-Kobayashi-Maskawa (CKM) 
phase P = a.Tg{—Vcdy*f,/VtdV^D [1]. Measurements of sin2/3 in tree-level b — >• sec transitions 
have established the first experimental evidence for CP violation in the B meson system [2,3]. 
The combined result of the latest measurements of sin2/3 in these transitions is in agreement 
with SM predictions [4]. 

The measurement of sm2j3 in decays dominated by a penguin loop-level b sqq transition 
is particularly interesting because of the sensitivity to physics beyond the Standard Model [5] . 
The B factory experiments have explored time-dependent CPV asymmetries in several such 
decays [6], namely 5° ^ [7, 8], 5° ^ r]'K° [7, 9], 5° ^ K+R-R^ [7, 10], B^ 

fRi [11] and 5° ^ R°y [12]. 

In this paper we present an update of our previous measurement of CP violation in 
— > RgTT^ decays, reported in reference [12]. The CKM and color suppression of the 
tree-level b — > suu transition leads to the expectation that this decay is dominated by a 
top-quark-mediated b — > sdd penguin amplitude, which carries a weak phase a.rg{VtbV*g) . If 
other contributions, such as the b — > suu tree amplitude, are ignored, the time-dcpcndcnt 
CPV asymmetry is governed by sin2/5 [13]. The bound on the deviation from sin2/3 due 
to Standard Model contributions with a different weak phase is ~ 0.2 from SU(3) flavor 
symmetry [14] and ~ 0.1 in model-dependent QCD calculations [15]. 

In addition to the CPV parameters we also present an update of our previous measure- 
ment of the branching fraction of 5° K^tt^ [16]. Existing experimental data on branching 
fractions for B Rn decays show a small discrepancy with respect to various calculations 
in the literature, the so-called 'Rtt puzzle' [17, 15]. In particular, the ratio of the B^ — > R^tt^ 
branching fraction to the other B — > Rtt branching fractions is about 2 standard deviations 
larger than inferred from isospin symmetry. Further experimental input might either resolve 
the puzzle or provide evidence for new physics [18]. 

2 THE BaBAR detector AND DATASET 

The results are based on a sample of 226.6 ± 2.5 million T{4S) BB decays collected in 
1999-2004 with the BABAR detector at the PEP-II e^e^ energy-asymmetric collider at the 
Stanford Linear Accelerator Center. The BABAR detector, fully described in [19], provides 
charged particle tracking through a combination of a five-layer double-sided silicon micro- 
strip detector (SVT) and a 40-layer central drift chamber (DCH), both operating in a 1.5 T 
magnetic field. Charged kaon and pion identification is achieved through measurements of 
particle energy- loss (dE/dx) in the tracking system and Cherenkov angle {6c) in a detector of 
internally reflected Cherenkov hght (DIRC). A segmented CsI(Tl) electromagnetic calorime- 
ter (EMC) provides photon detection and electron identification. Finally, the instrumented 
flux return (IFR) of the magnet allows discrimination of muons from pions. 
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3 ANALYSIS METHOD 



At the T{4S) resonance, time- dependent CPV asymmetries are measured by reconstructing 
the distribution of the difference of the proper decay times, At = tcp — ttag, where the tcp 
refers to the decay time of the signal 5° and ttag to the other B ( -Btag)- If and decay 
to a common CP eigenstate /, the At distribution follows 

Vgoi^t) = X [ 1 ± ( ,5/ sin (Am^At) - Cf cos (Am^At) ) ] , (1) 

where the upper (lower) sign corresponds to i?tag decaying as B^ (B^), t is the B^ lifetime 
averaged over the two mass eigenstates, Am^ is the mixing frequency, Cj is the magnitude 
of direct CP violation and S the magnitude of CP violation in the interference between 
mixing and decay. For the case of pure penguin dominance, we expect S^^^^o — sin2/3, and 

CrO'kO = 0. 

We search for B^ — > K^t:^ decays in BB candidate events selected using charged particle 
multiplicity and event topology [20]. We reconstruct K^g -k^tx" candidates from pairs 
of oppositely charged tracks. The two-track combinations must form a vertex with a 
consistency larger than 0.001, a tt+tt" invariant mass within 11.2 MeV/c^ (~ S.Scr) of the 
nominal mass [4] and a reconstructed decay length greater than five times its uncertainty. 
We form — > 77 candidates with an invariant mass 110 < m-y-y < 160 MeV/c^ from pairs of 
photon candidates in the EMC that are isolated from any charged tracks, carry a minimum 
energy of 50 MeV, and have the expected lateral shower shapes. B^ — > K^ti^ candidates 
are reconstructed from 7r° combinations and constrained to originate from the e^e~ 
interaction point using a geometric fit. We require that the consistency of the fit, which 
has one degree of freedom, be larger than 0.001. 

For each B candidate we compute two independent kinematic variables, namely the 
invariant mass tub and the missing mass m^iss — |?e+e- ~ where qe+e- is the four- 
momentum of the initial e^e~ system and qs is the four-momentum of the B^ — > K^^t:^ 
candidate after a mass constraint on the B^ is applied. Compared to the kinematic variables 

/S.E — E*Q — ^^/s and mES = \J\^ ~ P*b (where ^/s — \qe+e-\ and the asterisk denotes 
the center of momentum frame) that are traditionally used to select B decays in BABAR, 
the present combination of variables exhibits a smaller correlation and a better background 
suppression, in particular for modes with a relatively poor B energy resolution. We select 
candidates with ttib within 150 MeV/c^ of the nominal B^ mass [4] and 5.11 < mmiss < 
5.31 GeV/c^, which includes a sideband region for background characterization. 

To discriminate jet-like e'^e~ — > qq events (with q e {u, d, s, c}) from the more uniformly- 
distributed BB events we exploit the ratio L2/LQ of two Legendre moments defined as 
Lj = |p*|| cos^^*p, where p* is the momentum of particle i in the e+e" rest frame and 6* 
is the angle between p* and the thrust axis of the B^ candidate. We require L2/ Lq < 0.55, 
which suppresses the background by more than a factor 3 at the cost of approximately 10 % 
in signal efficiency. Finally, we require |cos^^| < 0.9, where 9*^ is the angle between the 
B^ candidate momentum and the momentum in the e^e~ rest frame. For B candidates 
the distribution of 0*^ follows P(cos^^) = 1 — cos^^^, whereas for continuum events it is 
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nearly flat. After all selections the average candidate multiplicity in events with at least 
one candidate is ~ 1.007. We select the candidate with the smallest on the 7r° mass as 
computed in the candidate vertex fit. 

For each K^n^ candidate we examine the remaining tracks and neutral candidates 

in the event to determine the decay vertex position and the flavor of i?tag- Using a neural 
network based on kinematic and particle identification information [21] each event is assigned 
to one of five mutually exclusive tagging categories, designed to combine fiavor tags with 
similar performance and At resolution. We parameterize the performance of this algorithm 
in a data sample (-Bflav) of fully reconstructed B^ D^*^^ / p'^ / af decays. The average 
effective tagging efficiency obtained from this sample is Q = Xlc^S'l-'-^^^^)^ ~ 0.288 ±0.005, 
where eg and w'^ are the efficiencies and mistag probabilities, respectively, for events tagged 
in category c e {1 ■ ■ ■ 5}. For the background the fraction of events (e^) and the asymmetry 
in the rate of B^ versus B^ tags in each tagging category are extracted from a fit to the data. 

To compute the proper time difference At the Stag vertex is inclusively reconstructed 
from the remaining charged particles in the event using the trajectory derived from the 
reconstruction of the Bcp candidate as a seed [20] . The time difference At and its uncertainty 
are extracted with a global fit to the T(45') — > B'^B^ decay tree that takes the information on 
the beam energy and the position of the interaction point (IP) into account. The position 
and size of the interaction region are determined on a run-by-run basis from the spatial 
distribution of vertices from two-track events. The uncertainty in the IP position follows 
from the size of the interaction region (about 200 fim horizontal and 4 /im vertical). 

Without additional constraints the single trajectory emerging from the S° — > ii'°7r° 
decay vertex provides insufficient information on the S° vertex position for a meaningful At 
measurement. To obtain the required resolution we constrain the sum of the two B lifetimes 
in the T{AS) B^B^ fit {tcp + ttag) to be equal to 2 r^o with an uncertainty \/2tbq- We 
have verified in a Monte Carlo simulation that this procedure leads to an unbiased estimate 
of At. 

For the ~ 40 % of events in which each of the two pion candidates from the K'^ decay does 
not have at least 4 SVT hits or for which a{At) > 2.5 ps or At > 20 ps, the At information 
is not used. However, since C can also be extracted from fiavor tagging information alone, 
these events still contribute to the measurement of C. 

We extract the signal yield and CPV parameters from an unbinned maximum-likelihood 
fit to mB, rrimiss, L2/LQ, cos^^. At and the fiavor tag variables. Exploiting sideband regions 
in data for the background and Monte Carlo simulation for the signal, we have verified that 
with the selection presented above these observables are sufficiently independent that we can 
construct the likelihood from the product of one dimensional probability density functions 
(PDFs). The PDFs for signal events are parameterized from either more copious fully- 
reconstructed B decays in data or from simulated events. For background PDFs we select 
the functional form from data in the sideband regions of the other observables, in which 
backgrounds dominate. We include these regions in the fitted sample and simultaneously 
extract the parameters of the background PDFs along with the CPV measurements. 

We obtain the PDF for the At of signal decays from the convolution of Eq. 1 with a 
resolution function TZ{5t = At — Attrue, c^At)- The resolution function is parameterized as 
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the sum of a 'core' and a 'tail' Gaussian, each with a width and mean proportional to the 
reconstructed UAt, and a third Gaussian centered at zero with a fixed width of 8ps [20]. 
We have verified in simulation that the parameters of TZ{St, a At) for B'^ Kgii'^ decays 
are similar to those obtained from the -Bfiav sample, even though the distributions of a^t 
differ considerably. Wc therefore extract these parameters from a fit to the -Bflav sample. We 
find that the At distribution of background candidates is well described by a delta function 
convolved with a resolution function with the same functional form as used for signal events. 
The parameters of the background function are determined in the fit. 

To extract the yield and the CPV asymmetries we maximize the logarithm of the extended 
likelihood 

(Ns+Nb) 

CiSf, Cf, Ns, Ns, fs, fB, a) = (^^^^ X 

n Wsfse^Vsixi, Sf, Cf) + NBfBe'hrBixuyi; d)] x 

n [^sil - fs)e'sVUxi;Cf)+NB{l - fB)e'BV'Bixf,d)] , 
ieii 

where / {II) is the subset of events with (without) At information. The probabilities Vs and Vb are 
products of PDFs for signal (S) and background (B) hypotheses evaluated for the measurements 
Xi = {ms, mjniss, ivi2/Lio, cos 0^, tag, tagging category} and yi = {At, cJAt}. Along with the CPV 
asymmetries Sj and Cf, the fit extracts the yields Ns and Nb, the fractions of events with At 
information fs and /b, and the remaining parameters, collectively denoted by a. These include all 
parameters of background PDFs and some parameters of the signal PDFs, such as the mean values 

of UlB and JTimiss- 



4 RESULTS 

Fitting the data sample of 9726 Kgir'^ candidates, we find Ns = 300 it 23 signal decays with 

5^o^o = 0.35 toil (stat) ± 0.04 (syst) 

and 

C^o^o = 0.06 ± 0.18 (stat) ± 0.06 (syst) . 

The total detection efficiency for K'^'^^ decays with K'^, 7r~^7T~ and vr'^ 77 is 0.34ib0.02 %. 

With the Kg and 7r° branching fractions taken from [4] and assuming equal production of charged 
and neutral B mesons at the T{4S) resonance, we obtain a branching fraction 

B{B° K^TT^) = ( 11.4 ± 0.9 (stat) ± 0.6 (syst) ) x 10"^ . 

The evaluation of the systematic uncertainties is described below. 

Figure 1 shows the so called s-Plot distributions [22] of the reconstructed mass and missing mass 
for signal B^ K'^t:^ candidates. The s-Plots are a statistical tool to extract the distribution 
of a particular discriminating variable by weighting events with their signal likelihood based on 
other discriminating variables. For comparison, figure 2 shows the unweighted distributions for a 
subsample enhanced in signal purity by selecting on L2/ and ?7imiss • 

The dashed and solid curves 
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Figure 1: Signal s-Plots for reconstructed mass (left) and missing mass (right). The curve 
represents the PDF used for signal events in the maximum likelihood fit. 




5.15 5.2 5.25 



5.35 5.4 
m...(GeV) 




Figure 2: Distributions for the reconstructed mass of the B candidate (left) and for the 
missing mass (right). To enhance the sample in signal purity we required L2/ Lq < 0.4 which 
reduces the efficiency for signal decays by ~ 25 % with respect to the selection described 
the text. For the distribution we required in addition m^iss > 5.25 GeV/c^. The 



m 



dashed and solid curves represent the background and signal-plus-background contributions, 
respectively, as obtained from the maximum likelihood fit. 
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Figure 3: s-Plot distributions of At for signal events with Btag tagged as (top) or 
(center), and of the asymmetry ^j^o^o(At) (bottom). The curves represent the PDFs for 
signal decays in the likelihood fit. 
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Figure 4: Contours for constant likelihood corresponding to a change in the likelihood with 
1, 2 and 3 units with respect to the minimum likelihood. The enclosed regions correspond 
roughly to 39%, 86% and 99% confidence levels. The star represents the Standard Model 
prediction. The large circle represents the boundary of the physically allowed region. 

indicate background and signal-plus-background contributions, respectively, as obtained from the 
fit, but corrected for the selection. 

Figure 3 shows the s-Plot distributions of At for B^- and S'^-tagged events, and of the asym- 
metry Aj^o^o(At) = [N^o — Ngo] I [-^so + -/V50] as a function of At. Figure 4 shows the contours 
for constant likelihood in the S — C plane. 

5 VALIDATIONS AND SYSTEMATIC UNCERTAIN- 
TIES 

The extraction of At with the IP-constrained geometric decay tree fit has been extensively tested 
on large samples of simulated K^tt^ decays with different values of C and S. To evaluate 

uncertainties due to the use of a resolution function extracted from the -Bflav sample we fit the 
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Kg-K^ samples with a resolution function extracted from simulated — > J/tp events. We 
assign a systematic uncertainty of 0.02 on Sj^o^^o and 0.01 on C^^o^o due to the use of the Bflav 
resolution function and other effects related to the /Si, reconstruction. We evaluate the effect of 
a possible misalignment of the SVT by introducing misalignments in the simulation and obtain a 
systematic uncertainty of 0.03 on Sj^o^o and 0.01 on C^o^o- We also consider large variations of 
the position and size of the interaction region, which we find to have negligible impact. We include 
a systematic uncertainty of 0.02 on both Sj^-o^o and C^^o^o to account for imperfect knowledge of 
the PDFs used in the fit. We assign a systematic uncertainty of 0.05 to C^^o^o due to possible 

asymmetries in the rate of B^ versus B^ tags in background events. 

In order to exclude any bias in the IP constraint At reconstruction that would be specific to 
data only, we examine a sample of approximately 1900 B^ — > J/tpK^ decays with J/^ — > 
J/^p e~^e~. In these events we determine At in two ways: by fully reconstructing the B^ decay 
vertex using the trajectories of charged daughters of the J/ip and the Kg mesons, or by neglecting the 
J/ijj contribution to the decay vertex and using the IP constraint and the Kg trajectory only. This 
study shows that within statistical uncertainties the IP-constrained At measurement is unbiased 
with respect to the more established technique, and that the values of Sj^^^j^o and Cjy^^o so 
obtained are consistent. In addition, we examine the pull of the At difference, assuming that the 
At uncertainties are fully correlated. We find that the pull distribution in the data is approximately 
10 % wider than that in the simulation. We therefore include an additional systematic uncertainty 
of 0.014 in Sj^o^o. Finally, we measure the B^ lifetime in B'^ K^tt^ decays and in IP-constrained 
B^ J/ip Kg decays and find that both agree with the world average. 

The detection efficiency for signal events is calculated with a Monte Carlo simulation based on 
the Pythia event generator [23] and GEANT4 detector simulation [24]. The efficiency of the Kg 
selection is calibrated with a large sample of inclusive Kg vr^vr^ decays. The tt" 77 efficiency 
is calibrated with e~^e~ t~^t~ events with r~ P~i't- The systematic uncertainty associated 
with the efficiency is 2.6% for K^ and 3.0% for 7r°. We assign additional systematic uncertainties 
of 1.2% for the L2/L0 cut, 2.0% for the selection on and a total of 2.0% for uncertainties in 
the signal PDFs. Finally, we include an uncertainty of 1.4% to account for unknown contributions 
from other BB decays and an uncertainty of 0.6% in the total number of T{AS) —>■ BB decays. 

6 SUMMARY 

In this paper we have reported preliminary results from a measurement of the branching fraction 
and time-dependent CPV asymmetries of B^ Kgir^ decays. The measured values of Sj^o^o and 
C^o^o are consistent with the Standard Model predictions. The measured branching fraction is 
consistent with measurements from other experiments [25]. 
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